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Abstract 

We consider the uplink of a multicell multiuser multiple-input multiple-output system where the 
channel experiences both small and large-scale fading. The data detection is done by using the linear 
zero-forcing technique, assuming the base station (BS) has perfect channel state information. We derive 
exact closed-form expressions for the uplink rate, symbol error rate, and outage probability per user, 
as well as a lower bound on the achievable rate. This bound is very tight and becomes exact in the 
large-antenna limit. We further study the asymptotic performance of the system in the regimes of high 
signal-to-noise ratio (SNR), large number of antennas, and large number of users per cell. We show that, 
at high SNRs, the system is interference-limited and hence, we cannot improve the system performance 
by increasing the transmit power of each user. Instead, by increasing the number of BS antennas, the 
effects of interference and noise can be reduced, thereby improving the system performance. We further 
demonstrate that, with very large antenna arrays at the BS, the transmit power of each user can be made 
inversely proportional to the number of BS antennas while maintaining a desired quality-of-service. 
Numerical and simulation results are presented to verify our analysis. 
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I. Introduction 

Multiple-input multiple-output (MIMO) technology can provide a remarkable increase in data 
rate and reliability compared to single-antenna systems. Recently, multiuser MIMO (MU-MIMO), 
where the base stations (BSs) are equipped with multiple antennas and communicate with 
several co-channel users, has gained much attention and is now being introduced in several 
new generation wireless standards (e.g., LTE-Advanced, 802.16m) [|2l. 

MU-MIMO systems have been studied from many perspectives including communication, 
signalling, and information theory in both downlink and uplink scenarios O, flU. For example, 
the downlink capacity of a single user in a MU-MIMO system has been investigated for optimum 
receivers in [5]. However, optimum receivers induce a significant complexity burden on the 
system implementation, especially for large MIMO configurations. Therefore, linear receivers, 
in particular Zero-Forcing (ZF), are of particular interest as low-complexity alternatives jH- 
[[8l. Note that all the above mentioned works have only investigated a single-cell scenario, 
where the effect of intercell interference has been neglected. However, co-channel interference, 
appearing due to frequency-reuse, represents an important impairment in cellular systems and can 
significantly degrade the system performance. Recently, there has been an increase of research 
work on the performance of MU-MIMO in interference-limited multi-cell environments ll9l- [[T3l . 
In fact, it has been shown that the capacity of the MIMO downlink can be significantly reduced 
due to intercell interference flU. 

Many interference cancellation or mitigation techniques have been proposed for multicell MU- 
MIMO systems, such as maximum likelihood multiuser detection ifTTI . BS cooperation lfT4l . and 
interference alignment [[T5l . These techniques, however, have a high implementation complexity. 
Very recently, there has been a great deal of interest in multicell MU-MIMO where the BSs are 
equipped with very large antenna arrays lfT6l - [|20l . With very large antenna arrays, the intercell 
interference can be significantly reduced with simple linear detectors, assuming the channel 
vectors between the BS array and the users are orthogonal enough (see e.g., [[T6l . [[T9l for a 
more detailed discussion of this assumption). In this context, approximations of the signal-to- 
interference ratio (SIR) when the number of BS antennas is very large are derived in [|T7l| . 
[|20l . To the best of the authors' knowledge, there is no analytical framework for analyzing the 
performance of uplink MU-MIMO systems in an interference-limited environment with a finite 
number of antennas. 
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Inspired by the above discussion, we hereafter analyze the performance of multicell MU- 
MIMO systems where a plethora of users simultaneously transmit data to a BS. The BS uses 
ZF to detect the transmitted signals. The paper makes the following specific contributions: 

• We derive exact closed-form expressions for the ergodic data rate, symbol error rate (SER), 
and outage probability of the uplink channel. We also derive a tractable lower bound on the 
achievable rate. These results enable us to explicitly study the effects of transmit power, 
intercell interference, and number of BS antennas on the system performance. 

• Next, we focus on the ZF receiver's asymptotic performance when the BS deploys a large 
antenna array. When the number of users per cell is fixed and the number of BS antennas 
grows without bound, intercell interference and noise are averaged out. However, when 
fixing the ratio between the number of BS antennas and the number of users, the intercell 
interference does not vanish when the number of antennas grows without bound. Yet, in 
both cases by using very large antenna arrays, the transmit power of each user can be made 
inversely proportional to the number of antennas with no reduction in the performance. 

Notation: The superscript H stands for conjugate transpose. [A] denotes the (i, j)th entry of a 
matrix A, and 7„ is the n x n identity matrix. The expectation operation and the Euclidean norm 
are denoted by E {■} and || ■ ||, respectively. The notation ^ means almost sure convergence. We 
use a = 6 to imply a and h have the same distribution. Finally, we use z ~ CN (0, E) to denote 
a circularly symmetric complex Gaussian vector z with zero-mean and covariance matrix E. 



In the following, we consider a multicell MU-MIMO system with L cells. Each cell includes 
one BS equipped with N antennas, and K single-antenna users (A^ > K). We consider uplink 
transmission, and assume that the L BSs share the same frequency band. Conventionally, the 
communication between the BS and the users is performed in separate time-frequency resources. 
However, this approach inherently reduces the spectral efficiency, and it is therefore better if 
several users communicate with the BS in the same time-frequency resource lfT3l . [|T7l . We 
assume that all users simultaneously transmit data streams to their BSs. Therefore, the x 1 
received vector at the /th BS is given by 



II. Multicell MU-MIMO System 



L 




(1) 
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where Gu G C 



NxK 



is the channel matrix between the /th BS and the K users in the ith cell, 



l-S-j Qlimk 



is the channel coefficient between the mth antenna of the Ith BS and the 



kth user in the zth cell; ^/p^Xi G C^^^ is the transmitted vector of K users in the ith cell (the 
average power transmitted by each user is pu); and Ui G C^^^ is an additive white Gaussian 
noise (AWGN) vector, rii ~ CJ\f {0,Im)- Note that, since the noise power is assumed to be 1, 
Pu can be considered as the normalized "transmit" SNR and hence, it is dimensionless. 

The channel matrix, Gu, models independent fast fading, path-loss attenuation, and log-normal 
shadow fading. Hence, its elements gumk are given by 



Qlimk 



limk ^ 



lik 1 



m 



1,2,. ..,iV 



(2) 



where humk is the fast fading coefficient from the A;th user in the ith cell to the mth antenna of 
the Ith BS. The coefficient humk is assumed to be complex Gaussian distributed with zero-mean 



and unit variance. Moreover, vAifc represents the path-loss attenuation and shadow fading which 
are assumed to be independent over m and to be constant over many coherent intervals. 

We assume that the BS has perfect channel state information (CSI). We further assume that 
the transmitted signals from the K users in the Ith cell are detected using a ZF receiver. As 
such, the received vector yi is processed by multiplying it with the pseudo-inverse of Gu as: 



ri = GjiPi = y/p^xi + ^^G\jGuXi + G^rii 



(3) 



where G\ = {G^Gu) G^ . Therefore, the kth element of ri is given by 



GiiXi + 



J k 



ni 



(4) 



ri,k = y/Kxi,k + [Gil 

where Xi^k is the kth element of Xi, which is the transmitted signal from the kth user in the Ith 
cell, while [A]f, denotes the A;th row of matrix A. From ([4]), the signal-to-noise-plus-interference 
ratio (SINR) of the uplink transmission from the kth user in the Ith cell to its BS is defined as 

Ik = ; — 2 ; — ; — 2- (5) 



J k 



Proposition 1: The SINR of the uplink transmission from the kth user in the Ith cell to its 
BS can be represented as 

d PuXk 



Ik 



PuZl + 1 



(6) 
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where and Zi are independent random variables (RVs) whose probability density functions 
(PDFs) are respectively given by 



X 



where Ai G 



PZi (^) = XI XI ^^^^ 

m=l n=l 



N-K 



X > 



— ^ z e^^' 



(n- 1)! 



z > 



(7) 
(8) 



D 



11 







-1) 



D 



1(1+1) 



D 



IL 



and Q {Ai) is the number of distinct diagonal elements of Af, yU./,2, lJ^i,e{Ai) the distinct 
diagonal elements in decreasing order; (^/) is the multiplicity of /i;,™; and Xm^n [Ai] is the 
(m, n)th characteristic coefficient of Ai which is defined in [[2T1 Definition 4]. 



Proof: Dividing the denominator and nominator of (|5]) by 

-2 

"J k 



we obtain 



Pn 



Ik 



^11 



where Yi 



~ IWgU.W 



. Since 



^ IIV.||2 



1 



-1 



fcfc 



(9) 



has an Erlang distribution 



with shape parameter N — K + 1 and scale parameter [l22l . Then, 



-2 



(10) 



Conditioned on 



Gi 



, Fj is a zero-mean complex Gaussian vector with covariance matrix Di 



which is independent of 



Gj, 



. Therefore, Yi ~ CM {Q.Du), where Du is a K x K diagonal 



matrix whose elements are given by [-D«j]^^. = Aifc- Then, II^jII^ is the sum of K {L — 1) 
statistically independent but not necessarily identically distributed exponential RVs. Thus, from 
||23l Theorem 2], we have that 

From (l9l)-(fTT]). we can obtain Q. 



Z,. 



(11) 
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III. FiNiTE-A^ Analysis 

In this section, we present exact analytical expressions for the ergodic uplink rate, SER, and 
outage probability of the system described in Section |Ill 



A. Uplink Rate Analysis 

From Proposition [H the uplink ergodic rate from the fcth user in the /th cell to its BS (in 
bits/s/Hz) is given by 



(i?,) = Ex„z, log2 1 + 



PnZl + 1 

^ ^ ^m,n Hi j^ 10g2 e 



cx> poo 



In 1 + ^ 

JO V PnZ + 1 



{ ^ (n-l)!(iV-AO<.-''^Vo Jo -\-puZ + l 



In (l + ^^lx^-^^e^^^-^e'^'^^dxdz. (12) 



We first evaluate the integral over x. By using [|24l Eq. (4.337.5)], we obtain 



(>{Ai)t„,{Ai) N-K y / A\ -n 1 

(^fc) = 



m=l 71=1 p=0 



(n-1)! (iV-if-|))! 7o 

q=l 



- {p^z + 1 



.N-K-p 



N-K-p 



PuZ + 1 
l^llkPn 



PukPn 

N-K-p-q 



e/'iifepuEi 



f^llkPu 

^n-lg-^M,„^^ (13) 



where Ei(-) is the exponential integral function ll24l Eq. (8.211.1)]. 

Theorem 1: The uplink ergodic rate from the kth user in the /th cell to its BS is given by 



g{Ai) t^{A,)N-K 
m=l n=l p=0 



'^m,n(A)/ir^(-l) 



N-K-p r- 



{n-l)\ {N-K-p)\ 



—el^UkP^Xn-l 



,N-K-'p 



1 



+ E IX^A r (n) [/ ( n, n + iV + 1 - - p - g, 



Pllk 

1 



(14) 



where U (■, ■, ■) is the confluent hypergeometric function of the second kind ll24l Eq. (9.210.2)], 



X^,„ (a, 6, a) = ^ ( J (-6) 



1=0 



q=0 



^q+l^rn-q 



QTi+i+l Qin-n-i 



_| i 1 

" ^ ^ a^a™-? (a/a + 1)^+' 



(15) 



Proof: See Appendix lAl ■ 
In practice, users are located randomly within cells, hence the large-scale fading factors for 
different users are different. This implies that all diagonal elements of Ai are distinct. Therefore, 
we provide the following corollary corresponding to this interesting case. 
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Corollary 1: If all diagonal elements of Ai are distinct, the ergodic rate in (fT4l) reduces to 

K{L-l]N^K y-fK(L-l) r \-l r / 1 1 1 1 \ 

/D \ 1 lln-1 [l^l.m- l^l,n) w^n- I ^ ^ ^ ^\ 



m=l p=0 



(A^-K-p)! (-1)- 



5=1 yilk 

Proof: For this case, substituting g (Ai) = K [L — 1), (^/) = 1, and 

K{L-l) , . _i 



(16) 



nil V ^/,™7 



into (fT4l) . and using the identity \J (1, a, x) = e^x^~'^r (a — 1, x), we can obtain (fT6l) . ■ 
In addition to the exact result given by Theorem [H we now derive an analytical lower bound 
on the ergodic achievable rate which is easier to evaluate. From (fT2)) . the uplink ergodic rate 
from the /cth user in the /th cell to its BS can be expressed as: 

= Ex„z, |log2 f 1 + Pu exp fin { ^\ \ \\\ (17) 



>log,(l+..exp(E.,,..{ln(-|^)})) (18) 

= log2 (1 + exp (Ex, {In (X^)} - Ez, {In (p,Z, + 1)})) (19) 

where we have exploited the fact that log2(l + aexp(x)) is convex in a; for a > along with 
Jensen's inequality. Therefore, we have the following new result: 

Proposition 2: The uplink ergodic rate from the /cth user in the /th cell to its BS is lower 
bounded by 

{Rk) > log2 (^1+Puf3ukexp (^ij{N -K+1) 

g(Ai)T^iAi) \\ 

) nui^rnsF'i in + 1,1, 1;2; -p^lii^rn) \ (20) 

m=l n=l / / 

where ^(x) is the Euler's digamma function [l24l Eq. (8.360.1)], and pFg{-) represents the 
generalized hypergeometric function with p,q non-negative integers ll24l Eq. (9.14.1)]. 

Proof: See Appendix |Bl ■ 
Remark 1: From we have that 

lim 7fc = ^^"^^^^„2 - (21) 
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The above result explicitly demonstrates that the SINR converges to a constant value when 
goes to infinity. This means that at high SNRs, we cannot improve the system performance by 
increasing the transmitted power of each user. The reason is that, when increases, both the 
desired signal power and the interference power increase. 

B. SER Analysis 

In this section, we analyze the SER performance of the uplink for each user. Let Ai^,. (s) be 
the moment generating function (MGF) of 7^. Then, using the well-known MGF-based approach 
II25I . we can deduce the average SER of M-ary phase-shift keying (M-PSK) as follows: 

Theorem 2: The average SER of the uplink from the kth user in the /th cell to its BS for 
M-PSK is given by 

^ ' 5'mpsk 



SERfc = - / M.,A^]d9 (22) 
where 6 = vr — -p, g^psK = sin^ (tt/M), and 

Proof: See Appendix O ■ 
It is also interesting to investigate the SER at high SNRs in order to obtain the diversity gain 

of the system performance. For this case (p^ — 00), by ignoring l/p^ in (|23l) . we obtain the 
asymptotic SER at high SNRs as 

SER- = ^ rM-(^)d9 (24) 



TT Jq Vsin 6 

where 

M-{s)=J2Y: E ''+'U.,™(^0(-l)^i^ok^;-;-^ . (25) 

m=l n=l p=0 ^ ^ ^ ^ P^^^-^ / 

This implies that at high SNRs, the SER converges to a deterministic constant, and hence the 
diversity order, which is defined as limp^^oo equal to zero. This phenomenon occurs 

due to the presence of interference. 

From (I22I) . we can see that to compute the SER we have to perform an integration over 6. 
To avoid this integration, we can apply the accurate approximation of ll26l on (l22l) . i.e.. 



. (I - 1) M,, (...) . \m.„ {^) . (I - 1) ) . (26, 

Clearly, the above expression is easier to evaluate compared to (l22l) . 
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C. Outage Probability Analysis 

The main goal of this section is to analytically assess the outage probability of multicell 
MU-MIMO systems. We recall that for the case of non-ergodic channels (e.g. quasi-static or 
block-fading), it is more appropriate to resort to the notion of outage probability to characterize 
the system performance. The outage probability, Pout, is defined as the probability that the 
instantaneous SINR, 7^, falls below a given threshold value 7th, i.e., 

Pout = Pr (7fc < 7th) . (27) 

With this definition in hand, we can present the following novel result: 
Theorem 3: The outage probability of transmission from the A;th user in the Zth cell to its BS 
is given by 

7th 



Pout = 1 - exp 



Ilk 



h h h h ^ - (V... + 7th/A..)"- • ^''^ 

Proof: See Appendix iDl ■ 
We recall that we are typically interested in small outage probabilities (e.g., 0.01, 0.001, ...). 
In this light, when 7th 0, we can obtain the following novel result: 

(iiAi)T„,{Ai) N-K ( Jl:k.Y -n , , 

poo _^ \^ \^ \PukJ_ , . . l^l,m ^ {n + p) 

m=\ n=l p=0 P ^ l^l,m + ltYvl Pllk) 

The above result is obtained by keeping the dominant term p = q and letting 7th 0. 

IV. Asymptotic {N 00) Analysis 

As discussed in Remark [TJ we cannot improve the multicell MU-MIMO system performance 
by simply increasing the transmit power. However, we can improve the system performance 
by using a large number of BS antennas. Due to the array gain and diversity effects, when 
increases, the received powers of both the desired and the interference signals increase. However, 
based on the orthogonal property of the channel vectors between the users and the BS, when N 
is large, the interference part can be significantly reduced even with a simple ZF receiver [[TtII . 
lfT9l . In this section, we analyze the asymptotic performance for large A^. 
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1) Fixed Pu, K, and N ^ oo: u 

Intuitively, when the number of BS antennas N grows large, the random vectors between the 
BS and the users as well as the noise vector at the BS become pairwisely orthogonal and hence, 
the interference from other users can be cancelled out. At the same time, due to the array gain 
effect, thermal noise cancels out too. This intuition is confirmed by the following analysis. 

Since Xk has an Erlang distribution with shape parameter N — K + 1 and scale parameter 
/3iik, Xk can be represented as 



where Zi, Z2, Ztv-a+i are independent, standard normal RVs. Substituting (l30l) into (|9l), and 
dividing the denominator and the numerator of 7^ by 2 (A^ — A' + 1), we obtain 



where (|3T| ) is obtained by using the law of large numbers, i.e., the numerator converges to 
Pu/3««fc/2, while the denominator converges to 0. The above result reveals that when the number 
of BS antennas goes to infinity, the effects of interference and noise disappear. Therefore, by 
increasing A^, the SINR grows without limit. 

2) Fixed pu, k = N/K, and N — )• oo.' This is an interesting asymptotic scenario since in 
practice, the number of BS antennas, A^, is large but may not be much greater than the number 
of users K. For this case, the property that the channel vectors between users and the BS are 
pairwisely orthogonal when N ^ 00 cannot be applied. In other words, the property that Hj^Hu 
converges point- wisely to an "infinite-size identity matrix" does not hold [18]. Therefore, the 
intercell interference cannot be cancelled out. Since Yi ~ CJ\f (0,Dii), it can be represented as 




2{N-K+1) 



(30) 



i=l 




— 00, as A^ — )■ 00 



(31) 



(32) 




(33) 



'We assume that when A'^ increases, the elements of the channel matrix are still independent. 
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By dividing the numerator and denominator of 7^ in (|33] ) by 2 (A^ — + 1), and using the law 
of large numbers and the trace lemma from ifTSl Lemma 1], i.e.E 



1 1 

—wfDuWi - —TrDu ^'0, as i^' 00 
K K 

we obtain a deterministic approximation of 7^ as follows 



a.s. 



7fc L 1 — — ' as ^ 00, and N/K = k. (34) 



It is interesting to note that the SIR expression (1341) is independent of the transmit power, and 
increases monotonically with n. Therefore, for an arbitrarily small transmit power, the SIR (|34|) 
can be approached arbitrarily closely by using a sufficiently large number of antennas and users. 
The reason is that since the number of users K is large, the system is interference-limited, so 
if every user reduces its power by the same factor then the limiting SIR is unchanged. 

Furthermore, from (l34l) . when k — )■ 00 (this is equivalent to the case N ^ K), the SIR 
7fc — 7- 00, as — 7- 00, which is consistent with (|3T1) . 

3) Fixed Np^, N — )■ oo.' Let = E^/N, where -Eu is fixed. From (|3T1) . we have 

Then, again using the law of large numbers and the trace lemma, we obtain 

7fc ^' A/fc^u, as N ^00, and fixed K (36) 

7. A^f^(l-V^) as AT ^ 00, N/K = k. (37) 

These results show that by using a very large antenna array at the BS, we can cut the transmit 
power at each user proportionally to 1/A^ while maintaining a desired quality-of-service. This 
result was originally established in [19] for the case when A^ ^ A' ^ 1 whereas herein, we 
have generalized this result to the regime where A^ ^ 1. Again, we can see that, when k tends 
to infinity, the two asymptotic results (l36l) and (l37l) coincide. 

Remark 2: We can see from (l36l) that when A^ grows without bound and K is fixed, the effects 
of interference and fast fading disappear. The only remaining effect is noise. Let us define the 



^ We assume that Du has a uniformly bounded spectral norm over K. 
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"massive MIMO effect" as the case where the system is ultimately limited by noised From 
(l37l) . when N grows large while keeping a finite k, the system is still limited by interference 
from the other cells. This interference depends mainly on k (the degrees of freedom), and when 
K — > oo, we achieve the massive MIMO condition. Therefore, an interesting question is: How 
many degrees of freedom k are needed in order to make the interference small compared to the 
noise (to obtain the massive MIMO condition)? Mathematically, we seek to find k, that satisfies 

l0g2 I 1 + A^fc^u (1 - \ ^ ^ ^^^.^^^ ^ ^ ^3g^ 

V ^u/«:Er=w;:|TrA. + iy 

where Rk,oo = log2 (1 + (3ukEu) is the ultimate rate which corresponds to the regime where 
N ^ K 1. We more closely address this fundamental issue via simulations in Section IVl 

Remark 3: When ^ Ji ^ 1 and = E^/N, using the property ^(x) = ln(x) + 1/x + 
0(l/x^) for X — 7- oo, and observing that the second term of the exponential function approaches 
zero, we can simplify (|20l) to get (i?^) > log2 (1 + f3iikE^) , which coincides with (|36l ). This 
implies that the proposed lower bound becomes exact in the large-antenna regime. 

V. Numerical Results 

In this section, we provide some numerical results to verify our analysis. Firstly, we consider 
a simple scenario where the large-scale fading is fixed. This setting enables us to evaluate the 
correctness of our proposed analytical expressions, as well as study the fundamental effects 
of intercell interference, number of BS antennas, transmit power of each user on the system 
performance. We then consider a more practical scenario that incorporates small-scale fading 
and large-scale fading including path-loss, shadowing, and random user locations. 
A. Scenario I 

We consider a multicell MU-MIMO system with 4 cells sharing the same frequency band. 
In all examples, except Fig. IH we choose the number of users per cell to be K = 10. We 
assume that all direct gains are equal to 1 and all cross gains are equal to a, i.e., /3uk = 1, and 
Ajfc = cfjVj ^ I, k = 1,2, ...,K (a can be thought of as an interference factor). Furthermore, 
we define SNR = Pu- 



' The massive MIMO effect was originally defined in |20 ] to occur when the system performance is limited by pilot 
contamination owing to the use of non-orthogonal pilots in different cells for the uplink training phase. However, here we assume 
perfect CSI, and we are considering a particular operating condition where the transmit power is very small (pu ~ l/N). 
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Figure \T\ shows the uplink sum rate in the /th cell versus SNR, at cross gain a = 0.1 and 
for different numbers of BS antennas N = 10, 20, 40, 60, 80 and 100. The simulation curves 
are obtained by performing Monte-Carlo simulations using ©, while the analytical and bound 
curves are computed via (fT4l) and (|20l ). respectively. We can see that the simulated and analytical 
results match exactly. As expected, when N increases, the sum rate increases too. However, at 
high SNRs, the sum rate converges to a deterministic constant which coincides with our analysis 
(|2T)) . Furthermore, a larger value of N makes the bound tighter. This is due to the fact that when 

grows large, things that were random before become deterministic and hence, the inequality 
used in (fTSi) will hold with equality (see Remark [3]). Therefore, the bound can very efficiently 
approximate the rate when N is large. It can be seen from the figure that, even for moderate 
number of antennas (N ^ 20), the bound becomes almost exact across the entire SNR range. 

The effect of interference for different is shown in Fig. |2l Again, the simulated and analytical 
results match exactly, and the bound is very tight. Interestingly, its tightness does not depend on 
the interference level but on the number of BS antennas. We can see that when the cross gain 
increases (and hence, the interference increases), the sum rate decreases significantly. On the 
other hand, the effect of interference decreases when the number of BS antennas grows large. 
For example, at a = 0.1, the sum rates are 3.76, 38.35, and 73.20 for A^ = 10, 50, and 500, 
respectively, while at a = 0.5, the sum rates are respectively 0.93, 19.10, and 50.80 for A^ = 10, 
50, and 500. This means that when increasing the cross gain from 0.1 to 0.5, the sum rates are 
reduced by 75.27%, 50.20%, and 30.60% for A^ = 10, 50, and 500, respectively. 

The power efficiency of large array systems is investigated in Fig. |3l Figure [3] shows the 
uplink sum rate in the /th cell versus A^ at a = 0.1, 0.3, and 0.5 for the cases of = 10 and 
Pu = 10/A^. As expected, with = 10/A^, the sum rate converges to a constant value when A^ 
increases regardless of the effects of interference, and with = 10, the sum rate grows without 
bound (logarithmically fast with A^) when A^ increases (see (|3TI) and (l36l)). 

Figure |4] shows the required number of degrees of freedom k to achieve 80% (77 = 0.8) 
and 90% (77 = 0.9) of a given ultimate rate Rk,oo, for a = 0.1, and a = 0.5. We use (l38l) to 
determine k. We can see that k increases with Rk,oo- Therefore, for multicell systems, the BS 
can serve more users with low data rates. This is due to the fact that when Rk,oo increases, the 
transmit power increases and hence, the interference also increases. Then, we need more degrees 
of freedom to mitigate interference. For the same reason, we can see from the figure that when 
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the interference factor a increases, the required k increases as well. 

In Fig. [51 the analytical SER curves are compared with the outputs of a Monte-Carlo simulator 
for different N. Here, we choose 4-PSK and a = 0.1. The "Analytical (Exact)" curves are 
computed using Theorem |2l and the "Analytical (Approx)" curves are generated using (|26|) . In 
addition, the high-SNR curves, generated via (|24l) . are also overlaid. It can be easily observed 
that the analytical results coincide with the simulation results. This validates the correctness 
of the proposed analytical expression. Furthermore, we can see that the "Analytical (Approx)" 
curves are accurate in all cases. As in the analysis of the sum rate, when the SNR is moderately 
large, the SER decreases very slowly and approaches an error floor (the asymptotic SER) due to 
interference, when SNR increases. Yet, we can improve the system performance by increasing 
the number of BS antennas. The effects of using large antenna arrays on the SER can be further 
verified in Fig. [6] where the SER is plotted as a function of N for different cross gains and 
4-PSK, at SNR = 10 dB. We can see that the system performance improves significantly when 
we increase A^. 

B. Scenario II 

We consider a hexagonal cellular network where each cell has a radius (from center to vertex) 
of 1000 meters. In each cell, K = 10 users are located uniformly at random and we assume 
that no user is closer to the BS than rh = 100 meters. The large-scale fading is modeled via 
Aifc = ^lik/ {riik/rhY, where zm, represents a log-normal RV with standard deviation of 8 dB, 
riik is the distance between the A;th user in the ith cell to the /th BS, and u is the path loss 
exponent. We choose u = 3.8 for our simulations. Furthermore, we assume that the transmitted 
data is modulated using OFDM. Let and be the OFDM symbol duration and useful symbol 
duration, respectively. Then, we define the net uplink rate of the kth user in the /th cell as follows: 



where B is the total bandwidth, and r is the frequency-reuse factor. Note that (l39l) is obtained 
by using the result in Proposition [H For our simulations, we choose parameters that resemble 
those of the LTE standard [fTTl : Ts = 71.4|Usec, and = 66.7/isec. We further assume that the 
total bandwidth of the system is 20 MHz. 

Figure |7] shows the cumulative distribution of the net uplink rate per user for different 
frequency-reuse factors r = 1, 3, and 7, and different number of BS antennas = 20, 100. 




(39) 
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We can see that the number of BS antennas has a very strong impact on the performance 
enhancement. The probability that the net uplink rate is smaller than a given indicated rate 
significantly decreases when N increases. We consider the 95%-likely rates, i.e., the rate is 
greater than or equal to this indicated rate with probability 0.95. We can see that 95%-likely 
rates increase with N; for example, with frequency-reuse factor of 1, increasing the number of BS 
antennas from 20 to 100 yields a 8-fold improvement in the 95%-likely rate (from 0.170 Mbits/sec 
to 1.375 Mbits/sec). Furthermore, when N is large, the random channel becomes deterministic 
and hence, the probability that the uplink rate is around its mean becomes inherently higher. 

When comparing the effects of using frequency-reuse factors, we can see that, at high rate (and 
hence at high SNR), smaller reuse factors are preferable, and vice versa at low rate. Furthermore, 
we can observe that the gaps between performances of different reuse factors become larger 
when increases. This is due to the fact that, when is large, the intercell interference can be 
significantly reduced; as a consequence, the bandwidth used has a larger impact on the system 
performance. Table H] summarizes the 95%-likely net uplink rates as well as their mean values. 

VL Conclusion 

In this paper, we analyzed in detail the uplink performance of data transmission from K single- 
antenna users in one cell to its A^-antenna BS in the presence of interference from other cells. The 
BS uses ZF to detect the transmitted signals. We derived exact closed-form expressions for the 
most important figures of merit, namely the uplink rate, SER, and outage probability, assuming 
that the channel between the users and the BS is affected by Rayleigh fading, shadowing, and 
path loss. 

Theoretically, when increases we obtain array and diversity gains, which affect both the 
desired and interference signals. Hence, from this perspective the performance is not dramatically 
affected. However, since when the number of BS antennas is large, the channel vectors between 
the users and the BS are pairwisely orthogonal, the interference can be cancelled out with a 
simple linear ZF receiver. As a consequence, by using a large antenna array, the performance of 
the multicell system improves significantly. Furthermore, we investigated the achievable power 
efficiency when using large antenna arrays at the BSs. Large antenna arrays enable us to reduce 
the transmitted power of each user proportionally to 1/A^ with no performance degradation, 
provided that the BS has perfect CSI. 
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Appendix 

A. Proof of Theorem [7] 

We first present the following lemma which is used to derive the closed-form expression of 
the achievable ergodic rate in Section IIII-AI 

Lemma 1: Let m,n a^h and a G M^. Then 

poo 

I x'^(ax + 6)"e-"^Ei(-(ax + 6))rfx =X„,„(a,6,a) (40) 
Jo 



where Xm,n (a, a) is given by (iTSl) . 

Proof: By change of variables z = ax + b in (|40|) . we have 



/ x"" {ax + hf e-"^Ei (- [ax + h))dx = / [z - z^'e-^YK {-z) dz 

Jo '^"^ Jb 

n \ / 



i=0 

where 



POO 

Jp (a, /i) = / z^e-^'Ei {-z) dz, /i > 0. 

^ a 



(42) 



Using partial integration and the fact that dEi (z) /dz = e^/z, the integral (l42l) for p G N can be 
evaluated as 



Jp (a, /i) = Ei {-z 



+ - / e"^" [p^P-^Ei i-z) + z^^^e^^l 

f^Ja 

p-1 

/Cp («, /i) + ^Jp^i («, /i) = ^ ( j) ' ICp^g (a, /i) + ( j) Vo («, /i) (43) 



p-i 

jrp_i (a,/i) = 

.=0 

where the second equality follows from [[241 Eq. (2.321.2)], where 



while the last equality is obtained using recursion. Finally, by using [[241 Eqs. (5.231.2) and 

(6.224.1)], Jo (a,/i) in ([43l) can be evaluated as 

1 e~"'^ 
Jo («, ;u) = Ei (- + 1) a) + Ei (-a) . (45) 



From (|4T]), (|43l), (I44l), and dMl), we arrive at the desired resuk (|40] 
Using Lemma m and fTT^, Eq. (39)], we can easily obtain ([141) . 
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B. Proof of Proposition |2] 

We can now evaluate the expectations in (fT9l) and we begin with Ex^. {In(Xfc)}, which can 
be expressed as 

n-N+K-l poo 

Ex, {In (Xk)} = :g / ln(a;)e-^/^-a;^-^(ia; = ^(iV - + 1) + ln(/3„fc) (46) 

(TV - Aj! Jq 

where we have used ll24l Eq. (4.352.1)] to evaluate the corresponding integral. 
The second expectation in (fT9l) requires a different line of reasoning: 



e{Ai) Tm{Ai) 

Ez, {In {p^Zi + 1)} = ^ J2 ^rn,n (A 



m=l n=l 



(n-l)!7o 

The integral X admits the following manipulations 



In {p^z + 1) e^i'^dz 



(47) 



where G""'" 



X, 



'-^2,2 



ai,...,ap 





1,1 








1,0 



1,3 





1 -n, 1,1 






1,0 



(48) 



denotes the Meijer's-G function ['24', Eq. (9.301)], and we have ex- 
pressed the integrand ln(l + 02;) in terms of Meijer's-G function according to ll28l Eq. (8.4.6.5)]. 
The final expression stems from ll24l Eq. (7.813.1)]. In addition, we can simplify (1481) as follows 



7- _ „ ij"'+^n^'^ 





-n,0,0 






0,-1 



r(n + 1)3^1 (n + 1, 1, 1; 2; -p^fxi, 



(49) 



where we have used 11241 Eq. (9.31.5)] to obtain the first equality and used [|28l Eq. (8.4.51.1)] 
to obtain the second equality. Combining (|47l) with (|49l) and after some basic simplifications, 
we get 

g(Ai) Tm{Ai) 

E^^ {ln(puZz + 1)} ^ Xrn,n (A) nni^m3Fi{n + 1, 1,1; 2; -p^Hi^ra) ■ (50) 

m=l n=l 

Substituting (|46l) and (l50l) into (fT9l) . we conclude the proof. 



C. Proof of Theorem |2] 

The MGF of 7^ is given by 

noo 

M,, is) = E {e-^^ } = / Ex, {e-^'' } pz, {z) dz. 

Jo 

Using the PDF of Xk given by ([7]), we have that 



(51) 
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x^-^ exp I -X I — + 



dx 



z + + fims, 
where the last equality is obtained by using [[24l Eq. (3.326.2)]. 
Substituting (|52|) into (ISTI) . and using ([8]), we get 



(52) 



m=l n=l 



n - 1 ! 



2^ + 1/Pu 



EE E 

m=l n=l p=0 



V 



Z + 1/pu + (SllkS 



dz 



n-\)\ \/3iikS + l/pn 



X 



+ 1 



(53) 



'0 vl/Pu + AzfcS 

where the last equality is obtained by using the binomial expansion formula. To evaluate the 
integral in (l53l) . we first express (^jj^-^:^^ ^0 terms of Meijer's-G function with the 
help of (M Eq. (8.4.2.5)], and then using the identity (M Eq. (2.24.3.1)], to obtain 



M^, (s) 



E E E 

m=l n=l p=0 



N -K + 1 
P 



I3uks 



n - 1)! \(3iikS + 1/pu 



X 



^1,2 
<-^2,l 



T{p) 



1^1,: 



j3iikS + 



l~n,l— p 




Finally, using ll24l Eq. (8.4.51.1)], we arrive at the desired result (l23l) . 



(54) 



D. Proof of Theorem |i] 

From Proposition [T] and dTT] ). we have 



Pr 



<7- 



th 



We can now express the above probability in integral form as follows: 



/>oo 

^'out = / Pr [Xk < 7th(^^ + l/Pu) I 
Jo 



The cumulative density function (CDF) of Xk can be shown to be equal to 

Fx,(x) = Pr {Xk < x) = 1 - exp " 



1 / X 



(55) 



(56) 



(57) 
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where we have used the integral identity Il24l Eq. (3.351.1)] to evaluate the CDF. Combining 
(l56l) with (|57l) . we can rewrite Pout as follows: 

7th(2; + l/Pu)^ 1 [ 7th(2 + 1/Pu) 

p=0 



ft.. = /_ I 1 - exp ( j—\ Y. ( ft, ) I P^,(-)''- 



1 - exp ( ) XI / ( ) + 



p=0 



N-K { 3^ 



p=0 









Pllk J 





/ s e{Ai)T,^{Ai) N~K ( J£!lV -n 

V PuPllkJ 1 1 n 



in-l)\ 

m=l n=l p=0 



+ ^"-1 exp (-Z + ^\] dz. (58) 

Jo V \^^l,m PllkJ J 



Using the binomial expansion, and [|24l Eq. (3.326.2)], we arrive at the desired result (|28l) . 
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Fig. 1. Simulated uplink sum rate, analytical expression and lower bound versus SNR (L = 4, K = 10, and a — 0.1). 
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Fig. 2. 
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Simulated uplink sum rate, analytical expression and lower bound versus the cross gain a (L 
10 dB). 
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Fig. 4. Degrees of freedom k required to aciiieve rjRk,oo versus Rk,ao (L = 4, a = 0.1 and 0.5). 
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Fig. 5. Simulated average SER and analytical expression versus the SNR for 4-PSK (L = A, K = IQ and a — 0.1). 
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Net Uplink Rate per User (Mbits/sec) 



Fig. 7. Cumulative distribution of the net uplink rate per user for frequency reuse factors 1, 3, and 7 (A'^ = 20, 100, SNR = 10 
dB, (Tshadow = 8 dB, and v — 3.8). 
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TABLE I 

Uplink performance of ZF with frequency reuse factors 1, 3, and 7, FORpu = lOdB, = 20, 100, 

CTshadow = 8dB, AND V = 3.8. 



Frequency Reuse 
Factor 


0.95-likely Net Uplink Rate per 
User (Mbits/sec) 


Mean of the Net Uplink Rate per User 
(Mbits/sec) 




A^=20 


A^=100 


A^=20 


A^=100 


1 


0.170 


1.375 


20.295 


48.380 


3 


0.150 


1.205 


10.668 


21.945 


7 


0.145 


1.350 


6.731 


12.546 



